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Keck Spectroscopy of distant GOODS Spheroidal Galaxies: 
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(N 1 ABSTRACT 

We analyze the evolution of the Fundamental Plane for 141 field spheroidal 
galaxies in the redshift range 0.2<z<1.2, selected morphologically to a magni- 
tude limit F850LP=22.43 in the northern field of the Great Observatories Origin 
Survey. For massive galaxies we find that the bulk of the star formation was 
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CN . completed prior to z = 2. However, for the lower mass galaxies, the luminosity- 

^ . weighted ages are significantly younger. The differential change in mass-to-light 

ratio correlates closely with rest-frame color, consistent with recent star formation 
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and associated growth. Our data are consistent with mass rather than environ- 
ment governing the overall growth, contrary to the expectations of hierarchical 
assembly. We discuss how feedback, conduction, and galaxy interactions may 
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explain the downsizing trends seen within our large sample. 
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Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution 
galaxies: formation — galaxies: structure - cosmology: observations 



1. Introduction 

The question of the assembly history of elliptical and SO (hereafter spheroidal) galaxies 
has come into sharp focus in recent years with the suggested discovery of their possible 
progenitors at redshifts above 2 (e.g., Chapman et al. 2003). However, the topic remains 
controversial both from a theoretical and observational standpoint. 



department of Physics and Astronomy, University of California at Los Angeles, Los Angeles, CA 90095; 
ttreu@astro.ucla.edu; Hubble Fellow 

2 Caltech, Astronomy, 105-24, Pasadena, CA 91125, USA; rse@astro.caltech.edu, ting@caltech.edu 

3 Department of Astronomy, Yale University, New Haven, CT 06520; dokkum@astro.yale.edu 



- 2 - 



Semi-analytic models based on the current ACDM structure formation picture predict 
for spheroidals remarkably different assembly rates than observed (Benson et al. 2002), sug- 
gesting the effects of merging on stellar mass assembly are poorly understood. Consequently, 
the implication of recent claims for an abundant population of old massive systems at z ~2 
are unclear (Glazebrook et al. 2004; Cimatti et al. 2004). 

Exactly how the slow evolutionary trends observed in cluster spheroidals (Kelson et al. 
1997; Bender et al. 1998; Stanford, Eisenhardt & Dickinson 1998; van Dokkum & Stanford 
2003) contrast with those measured in the field (Treu et al 1999,2001b, 2002; van Dokkum 
& Ellis 2003; Gebhardt et al. 2003; Rusin et al. 2003; Bernardi et al. 2003; Yi et al. 2005; 
Jimenez et al. 2004) is also unclear. The environmental dependencies do not appear to be as 
large (van Dokkum et al. 2001) as predicted by some renditions of hierarchical models (e.g. 
Diaferio et al. 2001; de Lucia et al. 2004). And those trends indicative of early assembly 
appear to be strongly at variance with the rapid decline in the red luminosity density claimed 
from large photometric samples (Bell et al 2004). 

In this Letter we present the first results from a comprehensive Keck survey of 141 field 
spheroidals selected in the northern Great Observatories Origins Survey (GOODS) field. The 
spectroscopic signal to noise enables us to measure accurate stellar velocity dispersions for 
each spheroidal and, utilizing the Fundamental Plane (FP, Djorgovski & Davis 1987; Dressier 
et al. 1987), the dynamical mass and mass-to-light ratio. As these masses are closely linked 
to those of the underlying halos (Treu & Koopmans 2004), we establish a more complete 
picture of the assembly history than from studies based solely on photometrically-derived 
stellar masses. Our sample achieves a high degree of statistical completeness and is almost 
an order of magnitude larger than any previous equivalent study. This enables us to examine 
the rate of assembly as a function of dynamical mass and environmental density, both key 
variables for testing the currently-popular hierarchical picture for galaxy formation. 

We assume throughout a cosmological model with Qm = 0.3, fl\ = 0.7 and H =65 
km s" 1 Mpc" 1 . All magnitudes are in the AB system (Oke 1974), unless otherwise noted. 

2. Data 

In order to compare properly with local E+S0 galaxies, our parent sample was selected 
morphologically using visual classification to a F850LP magnitude limit of 22.43 (AB). In- 
ternal tests (Treu et al. 2005, in preparation, hereafter T05) show that our classification is 
self-consistent. The surface density of spheroidals is in agreement with that found by compa- 
rable studies (e.g. 1.32±0.09 arcmin -2 vs. 1.23±0.10 arcmin -2 measured by Im et al. 2002 
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to the same depth), showing that our classification is consistent with previous work and that 
the surface density is comparable with that expected from the local luminosity function. 

High signal-to-noise spectra were obtained during April 1-5 2003 using the DEIMOS 
spectrograph at the Keck-II Telescope. Conditions were excellent with seeing in the range 
0%-ff!8 throughout. The 1200 grating blazed at 7500A was centered at 8000 A. Total 
exposure times ranged between 14,400 and 38,900 seconds. The DEIMOS spectra were 
reduced using the IDL pipeline developed for the DEEP2 Survey (Davis et al 2003). 

Stellar velocity dispersion were derived for 141 spheroidals (out of 163 spectroscopically 
targeted) using the methods described in Treu et al. (1999,2001a). To estimate systematic 
errors, a subsample was analyzed independently using methods described in van Dokkum 
& Franx (1996). The independent determinations agree within 1% on average (12% rms), 
indicating that systematic differences due to methodology are negligible. Comparisons with 
previously-published dispersions in the Hubble Deep Field (van Dokkum & Ellis 2003) as 
well as repeat DEIMOS measures for targets on different masks show that our estimated 
errors (~ 10%) are reliable. Surface photometry through r 1 / 4 model fitting was derived 
from the images using the GALFIT package (Peng et al. 2002). 

A comparison of the subsample with measured velocity dispersions with the parent 
sample of morphologically selected spheroidals shows it to be representative with respect to 
color, luminosity and redshift distribution. We conclude that our sample is representative 
of morphologically selected spheroidals. The catalog, a description of the sample, related 
selection effects, observations and data reduction are given in T05. 

In order to test for environmental trends, local projected densities were derived for indi- 
vidual spheroidals using a modification of the algorithm described in Dressier (1980) and the 
GOODS-N photometric catalog of Bundy, Ellis & Conselice (2004; hereafter B04). The B04 
catalog includes spectroscopic redshifts from this survey and the Keck Team Redshift Survey 
(Wirth et al. 2004), with relatively uniform spatial sampling, supplemented by photometric 
redshifts (for the remaining 45% of the objects). For each spheroidal, we identify nearby 
objects brighter than My = —19.8, allowing for passive evolution. Then, moving away on 
the plane of the sky, we add the probability that each neighbor lies within 3000 km s _1 until 
a total number 5 is reached. This defines the local density S 5 . We chose 5 nearest neighbors 
(rather than the usual 10) to enhance sensitivity to more modest overdensities, resulting in 
a poisson noise component of ~0.2 dex. The partial use of photo-z brings the total error on 
E 5 to 0.3 dex per galaxy (as calculated with simulations) . 
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3. Results 
3.1. Evolution of the Fundamental Plane 

The Fundamental Plane (FP) is defined as: 

logi^e = « log a + P SB e + 7- (1) 

where Re and SBe are the effective radius and surface brightness and a is the central velocity 
dispersion. A physical interpretation of the evolution of the FP follows from defining an 
effective (dynamical) mass, 

M = —^. (2) 

Redshift variations of the slopes a, (5 and intercept 7 can be interpreted via evolution of the 
stellar populations. If a and R e and Eq. 2 do not evolve, for the ith galaxy: 

T^logi^-alog^-^SB*,, (3) 

and the offset from the FP ( A7* = 7* — 7) is related to the offset of the Mj L by 

The scatter of the FP (~ 0.1 in log_R e ; Bernardi et al. 2003) is thus a measure of the 
homogeneity of the stellar population (see Treu et al. 2001b for more discussion). 

Figure 1 shows the offset from the relation for the local Coma Cluster 1 using Eq. 4 for 
all 141 GOODS-N spheroids coded by dynamical mass (Eq 2). Evolutionary tracks for single 
burst stellar populations (Treu et al. 2001b) at various formation redshift (z { = 1,2,5) are 
shown for comparison. The most massive spheroidals uniformly follow a passive evolutionary 
trend consistent with an early epoch of formation, as found for massive cluster galaxies (e.g. 
van Dokkum & Stanford 2003). However, those with smaller masses are systematically 
younger (as hinted previously by smaller samples; Treu et al. 2002, van der Wei et al. 2004) 
and show increased scatter, implying a lower redshift of formation or secondary episodes of 
activity which rejuvenated a previously old population (e.g. Trager et al. 2000; Treu et al. 
2002; van Dokkum & Ellis 2003). 

We can verify that recent star formation is responsible for the trends in Figure 1 by 
comparing M/L with independent diagnostics of recent stellar formation, such as the rest- 
frame color (Figure 2). If we assume that the stellar M*/ ' L B can be obtained as AlogM/L B 



1 We use Coma as the local reference for both cluster and field to minimize systematic uncertainties related 
to filter transformations, distance determination and selection effects; sec discussion in Treu et al. 200f b. 
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+ logM*/LB, z =o with M*/L B ,z=o — 7.3M /L B ,0 (Treu & Koopmans 2004), we can reproduce 
the locus for an ageing single burst population with a Salpeter IMF (red line in Figure 1). 
The observed correlation is fairly close to this model, although the scatter is significant, 
presumably as a result of a mixture of primary and secondary episodes of star formation, 
and, possibly, of metallicities. The analysis of spectral diagnostics such as Balmer absorption 
lines (T05) is also consistent with the trend in M/L being due to recent starformation. 

3.2. The role of mass and environment 

Our large and homogeneous sample allows us to take a significant step forward with 
respect to previous studies by exploring further the trends of M/L with dynamical mass 
introduced in §3.1 and by examining the correlations between M/L and local galaxy density, 
£5. Dynamical mass can be an effective proxy for the total halo mass (Treu & Koopmans 
2004) and local density (§ 2) for environmental effects. Therefore, we can use these correla- 
tions to study the role of halo mass and environment in the context of contemporary models 
of structure formation. To this aim, we will consider the quantity 5A\ogM/L B , defined 
as the difference between Alog M/L B for an individual spheroidal and that expected for 
old massive cluster galaxies at that redshift (i.e. Alog M/L B = —0.46,2, van Dokkum & 
Stanford 2003). Negative values of 5 A log M/L B indicate a smaller M/L B than that of a 
massive cluster spheroid. 

Figure 3 shows SA log M/L B as a function of galaxy mass; a clear trend is evident with 
lower mass galaxies displaying younger luminosity-weighted ages. Although the two axes 
are not independent (a larger mass implies a larger M/L neglecting errors on L B ), both the 
small error bars (typically < 0.1 dex) and correlation with independent diagnostics (Figure 2; 
T05) indicate this trend does not arise from correlated uncertainties. Because of age and 
luminosity limits, a parent mass-selected sample must generally occupy a triangular region 
between the horizontal blue line and the lowest observed 5A\ogM/L B at any given mass. 
We thus conclude that the luminosity-weighted age of spheroidals increases with dynamical 
mass, a result we statistically confirm using detailed Monte Carlo simulations in T05. 

We illustrate the effects of secondary bursts of star formation in Figure 3 by considering 
a family of two burst models, arranged by adding a young stellar component ( aged 1 Gyr) 
to an earlier population (10 Gyr) equal for simplicity to the dynamical mass. For example 
a secondary burst of 1O 9 -1O 1O M would reproduce the observed trend. 

Figure 4 shows 5 A log M/L B as a function of local density £5. No statistically significant 
correlation is present, although the dynamic range in GOODS-N does not include massive 
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clusters. Within our larger error bars, our data are consistent with the trend observed 
locally in the Sloan Digital Sky Survey (SDSS; Bernardi et al. 2003), which in turn can 
be explained in terms of mass differences in the context of environment-dependent mass 
functions. More pronounced is the difference with massive galaxies in the cores of rich 
clusters (5 A log Mj Lb = 0), although the poor statistics available for clusters at comparable 
redshifts and the different selection criteria prevent a detailed comparison. 

4. Discussion 

We have shown that mass is a key variable determining the star formation history of 
field spheroidals. The most massive examples (> 10 11 ' 5 M ) formed most of their stars 
before z ~ 2, while less massive ones (< 10 11 M Q ) show greater scatter and systematically 
younger stellar ages. The secondary activity necessary to reproduce our observations (up to 
~ fO 10 M Q over the past few Gyr) could be explained by accretion of satellites or gas-rich 
mergers. Regardless of the process, substantial recent growth is implied in many cases for 
galaxies less massive than fO 11 M Q . 

Such "downsizing" (i.e. stellar growth that moves systematically to lower mass systems 
as the Universe expands; Cowie et al. f 996) appears contrary to the expectations of hier- 
archical models of galaxy formation (Diaferio et al. 2001). In this widely-discussed picture, 
massive spheroidals assembled recently via major mergers with associated star formation. 
Strong environmental trends are also expected in such a picture (e.g. De Lucia et al. 2004), as 
a result of the fact that evolution is accelerated in these "biased" regions. We do not find any 
environmental trend within the range of densities covered by the GOODS-N field, although 
a larger sample of data will be needed to confirm this and properly compare with dense 
clusters, by studying the bivariate distribution (mass and local density) of 5 A log Mj Lb- 
Out observations may be reconciled with the hierarchical paradigm via the introduction 
of more detailed physics that regulates star formation in collapsing systems. If, above a 
certain halo mass, feedback from Supernovae or AGN, or thermal conduction (Benson et al. 
2003; Nagashima et al. 2004; Granato et al. 2004), heat or remove gas and prevent it from 
cooling back onto the galaxy - thus inhibiting star formation - then the most massive galax- 
ies will appear old even if they assembled dynamically at recent times, e.g. via collisionless 
mergers. Alternatively, star formation could be anticipated in massive system, triggered 
by galaxy interactions (Menci et al. 2004). Although infalling satellites or minor mergers 
may be frequent, they are likely to influence only the lower mass systems, as observed. The 
extended timescale necessary for lower mass mergers to relax and present smooth morpho- 
logical systems may also play a key role (Hernandez & Lee 2004). Regardless of the detailed 
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explanation, it is clear that the mass-dependent growth of spheroidals is more complex than 
previously thought. It will be interesting to see if the numerical models with these additional 
mechanisms can reproduce quantitatively our result. 
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Fig. 1. — Top panels: edge-on view of the FP in redshift bins. The solid line represents the 
FP for the Coma cluster. Bottom panel: offset from the Coma FP, coded by dynamical mass 
demonstrating the stronger evolution for smaller masses. Typical errors are < 0.1 dex on 
the vertical axis. Dotted lines represent passive evolutionary trends for single burst models 
at various formation epochs (z{). 
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M./L B (M /L B0 ) 

Fig. 2. — Rest frame (B-V) color as a function of the stellar mass to light ratio, obtained 
assuming a local value of 7.3 M /Lb iQ . The red curves are the tracks followed by single 
burst stellar populations with Salpeter IMF (ages 0.5 - 12 Gyrs) and various metallicities. 
The size of symbols represents galaxy mass as in Figure 1. 
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Fig. 3. — M/Lb offset with respect to the redshift-dependent behavior of old massive cluster 
spheroidals, plotted against dynamical mass M. The blue horizontal line represents the locus 
seen for massive cluster galaxies. Red curves illustrate the effects of secondary activity (1 
Gyr ago) applied to a 10 Gyr old population of mass M. The blue (dotted) diagonal lines 
demonstrate the effect of the magnitude limit for each redshift bin. 
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Fig. 4. — Offset of field spheroidals from the cluster trend (5AM/ Lb) as a function of local 
projected density S 5 , coded by mass as in Figure 1 (The error on the local density is ~ 0.3 
dex). Solid squares represent the averages in each density bin; large error bars represent the 
scatter, small error bars the uncertainty in the mean. The black rectangle represents the 
value observed for massive cluster spheroidals, the diagonal dashed line the local trend from 
the SDSS Survey (Bernardi et al. 2003). 



